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' Abstract 
=^ ' 

A bicriteria approximation algorithm is presented for the unrooted travehng repairman prob- 
' lem, realizing increased profit in return for increased speedup of repairman motion. The algo- 

04 . rithm generalizes previous results from the case in which all time windows are the same length 

to the case in which their lengths can range between 1 and 2. This analysis can extend to any 
, range of time window lengths, following our earlier techniques This relationship between 

Q ' repairman profit and speedup is applicable over a range of values that is dependent on the cost of 

putting the input in an especially desirable form, involving what are called "trimmed windows." 
For time windows with lengths between 1 and 2, the range of values for speedup s for which our 
analysis holds is 1 < s < 6. In this range, we establish an approximation ratio that is constant 
for any specific value of s. 

> 
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m : 1 Introduction 

In this paper we present an approximation algorithm for a practical time-sensitive routing problem, 
the unrooted traveling repairman problem with time windows. The input to this problem is a speed 
at which a repairman can travel and a list of service requests. Each service request is located at a 
node in a weighted metric graph, whose edges give the travel distance between nodes. Each service 
^ ' request also has a specific time window during which it is valid for service. The goal of the problem 

' is to plan a route called a service run that, starting at any service request at any time, visits as 

many service requests as possible during their respective time windows. 

Because the problem is NP-hard, our only hope for an efficient approach seems to be an ap- 
proximation algorithm. In the real world, a repairman may have some flexibility in choosing speed. 
As a consequence, our earlier approximation algorithms [12] and this paper are parameterized by 
speedup s, so that we can characterize how much closer to optimal the repairman can do if he or 
she travels a factor of s faster than a hypothetical repairman traveling along an optimal route at 
the baseline speed. This type of approximation based on resource augmentation is well known in 
the scheduling community as shown by Bansal et al. [3], Kalyanasundaram and Pruhs [13j, and 
Phillips et al. [E]. 

The algorithms in this paper build on our earlier work [10\ [TT] , in which we introduced the first 
polynomial-time algorithms that give constant approximations to the traveling repairman problem 
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when all the time wmdows are the same length. As a counterpoint to the repairman problem, 
we also introduced the speeding deliveryman problem in [101 El], with an alternative optimization 
paradigm, namely speedup. The input to the speeding deliveryman problem is the same as the 
input to the traveling repairman problem, but the goal is to find the minimum speed necessary 
to visit all service requests during their time windows and thus collect all profit. In pJLj we also 
gave constant-factor polynomial-time approximation algorithms for both problems when the time 
windows have lengths in some fixed range. 

In both the repairman and deliveryman problems, our algorithms |10l fTT] rely on trimming 
windows so that the resulting time windows are pairwise either identical or non-overlapping . We 
trim time windows by repeatedly making divisions in time after a fixed amount of time has passed, 
starting at a specified time. We define a period to be the time interval that starts at a particular 
division and continues up to the next division. When time windows are unit length, we choose a 
period length of .5 time units. Because we define periods so that no window starts on a period 
boundary, each time window will completely overlap exactly one period and partially overlap its 
two neighboring periods. Trimming then removes those parts of each window that fall outside of 
the completely overlapped period. 

In this simpler case where time windows all have the same length, the penalty for trimming the 
repairman is a reduction by a factor of 1/3 in the number of requests serviced, and the penalty 
for the deliveryman is an increase by a factor of 4 in the speed needed to service all requests. In 
|12| . we showed that, for unit time windows, a spectrum of performance is possible between these 
two extremes. For some speedup s greater than 1 but less than 4, we showed how to achieve an 
increase in the number of serviced requests, proportional in some sense to s. The approximation is 
also a function of graph property 7, where 7 = 1 for a tree and 7 is no more than 2 + e for a metric 
graph, for any constant e > 0. A more complete explanation of 7 is given in Sect. [2j 

In this paper, we extend our algorithms and analysis to the more challenging case in which time 
windows have lengths in some fixed range, specifically between 1 and 2. We present an algorithm 
that finds approximations parameterized by speedup s and property 7. To prove these approxi- 
mation bounds, our analysis establishes and takes advantage of the existence of an ensemble of 
runs that move backward and forward along the path of an (unknown) optimal run, similar to our 
work in |12j . These runs are analyzed based on several different starting points for trimming. To 
handle windows of different lengths (i.e., between 1 and 2), we orchestrate several complementary 
trimming schemes, run our approximation algorithm on each combination, and choose the best 
result. 

On the surface, the approach we use to orchestrate trimming schemes is similar to our approach 
in [11], which extended our earlier approximation algorithms from [TO] to achieve a constant ap- 
proximation on time windows whose lengths were between 1 and 2. However, the similarity of the 
algorithms belies the fundamental difference in the analysis, whose complexity increases by at least 
an order of magnitude in the process of uniting speedup with non-uniform time windows. The 
key to our algorithm remains using a different period length for each trimming scheme, with each 
subsequent scheme using a progressively longer period length. Intuitively, by selecting the most 
profitable run found in any scheme, the algorithm adapts to different distributions of window sizes. 
If most of the windows are short, a scheme of trimming to shorter lengths will be effective. If 
most of the windows are long, a scheme of trimming to longer lengths will be effective. Because 
the output of each trimming scheme is a set of trimmed windows of equal length, our speedup 
algorithms from [12j can then be applied directly. As with the case of no speedup, we bound the 
approximation guarantee of our algorithm by accounting for a variety of distributions of windows, 
but the tool needed to bound each distribution is now a considerably richer set of hypothetical 
runs. 
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The major contributions of this paper are two additional techniques needed to extend the 
analysis for speedup on unit-time windows to windows with non-uniform length. The first technique 
is a significantly more complex design of ensembles to achieve good coverage, using a greater variety 
of runs, some of which have longer repeating patterns. Once we select an appropriate ensemble, 
we use a symbolic description of the coverage of the runs in the ensemble to demonstrate good 
coverage for all speedups in the range of speedups under consideration. The second technique is 
an approach for designing and coordinating together the different bounds of approximation as a 
function of speedup for different window lengths. Using averaging arguments, we will show that any 
convex combination of the approximation guarantees for each trimming scheme is a lower bound 
on the profit of the best run produced by our algorithm. For each range of speedups in question, 
we determine the best choices of weightings for a convex combination of the approximation bounds 
we have found. By using the best convex combination of bounds from each scheme, we guarantee a 
good bound of approximation. The details of these techniques are given for the case when window 
size is between 1 and 2, but other ranges of window size can be accommodated in a similar way. 

As a result, we can still produce polynomial-time approximation algorithms with constant- 
factor approximations for a given s over a significant speedup range. Our process of combining 
together different approximation bounds, as a function of the speedup s, gives a final result in 
Table [1] that is more involved than our results in [12] . The ratio has more piecewise ranges and its 
inverse is primarily nonlinear, even though the inverse of the ratio in each range is fairly close to a 
linear function. Note that approximation ratios are typically defined to be at least 1, and so these 
approximation ratios will give the reciprocal of the fraction of profit collected at a given speedup. 
For ease of presentation, most of the analysis in this paper will instead be in terms of the fraction 
of profit collected. 
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Table 1: Approximation ratios for speedup s when time window lengths are between 1 and 2. 

Our results are recent developments in time-sensitive routing problems, which have received a 
lot of attention from the algorithms community in the last decade. As with our particular problem, 
these problems typically identify the locations to be visited and the cost of traveling between them 
as the nodes and edges, respectively, of a weighted graph. For example, the orienteering problem 
considered by Arkin et al. [Ij, Bansal et al. [2], Blum et al. [5], Chekuri et al. p'J, and Chen and 
Har-Peled [9] seeks to find a path that visits as many nodes as possible before a global time deadline. 
The deadline traveling salesman problem which was also considered by Bansal et al. [2] generalizes 
this problem further by allowing each location to have its own deadline. Our traveling repairman 
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problem can be viewed as a further generalization from a deadline to a time window. 

A great deal of work by Bansal et al. |2], Bar- Yehuda et al. [3], Chekuri and Kumar [8], Karuno 
et al. 113], Tsitsiklis |17j . and the authors [TT] has been done on the traveling repairman problem, 
although much of the preceding literature, including that from Bansal et al. [2] and Bar- Yehuda 
et al. [3], considers the rooted version of the problem, in which the repairman starts at a specific 
location at a specific time. 

For general time windows in the rooted problem, an 0(log^ n)-approximation is given by Bansal 
et al. [2]. An 0(log L)-approximation is given by Chekuri and Korula for the case that all time 
window start and end times are integers, where L is the length of the longest time window. In 
contrast, a constant approximation is given by Chekuri and Kumar [8], but only when there are a 
constant number of different time windows. Our earlier work [11] and work by Chekuri and Korula 
give 0(log L')-approximations to the unrooted problem with general time windows, where D 
is the ratio of the length of largest time window to the length of the smallest. Polylogarithmic 
approximation algorithms for the directed traveling salesman problem with time windows have also 
been given by Chekuri et al. [7J and Nagarajan and Ravi [15] . 

2 Trimming Time Windows and the Associated Loss 

In our earlier work |lHll2j. we trimmed time windows of unit length by first making divisions in time 
every .5 time units, starting at time 0. We generalize the process for time windows with different 
lengths by instead making divisions every a time units, where a = .5, .75, or 1 for the range of time 
window lengths [1,2). Define a period to be the time interval that starts at a particular division 
and continues up to the next division. In the case of unit time windows, each time window will 
completely overlap exactly one period and partially overlap its two neighboring periods, because 
we allow no window to start on a period boundary. In the case of longer time windows, there 
will be different patterns of overlapping. If a window completely overlaps with only one period, 
trimming will remove those parts of each window that fall outside of the completely overlapped 
period. If a window completely overlaps with more than one period, one trimming scheme will 
remove all those parts of each such window that fall outside of the first completely overlapped 
period. Another separate trimming scheme will remove all those parts of each such window that 
fall outside of the second completely overlapped period. For long period sizes, some time windows 
may not completely overlap any full periods and will vanish in the process of trimming. In each 
case, because periods do not overlap and a time window is trimmed to at most one period, trimmed 
time windows will be pairwise either identical or non-overlapping. Our repairman algorithm in [11] 
identifies a variety of good paths inside each separate period and then uses dynamic programming 
to select and paste these paths together into a variety of longer good paths and ultimately a good 
service run for the whole problem. 

To describe our results for both trees and general metric graphs, we use the graph property 7, 
where 7 = 1 for a tree, derived in our earlier paper [TT], and 7 < 2 + e for a metric graph, derived 
by Chekuri et al. [7], for any constant e > 0. To describe the running time for these repairman 
algorithms we use r(n), where T{n) is O(n^) for a tree and 0(n'^(^/'^ for a metric graph. The value 
of 7 and running time of r(n) are dependent on the approximation bounds for finding maximum 
profit paths within a specific period on a specific class of graph. Although the available results only 
give 7 values and r(n) running times for trees and metric graphs, other classes of graphs, such as 
outerplanar or Euclidean graphs, may have intermediate values of 7 and r(n). 

In [TT] we showed that, for unit time windows with no speedup, the reduction due to trimming 
still allows us to visit at least 1/(87) of the pre-trimming optimal and, with a speedup of 4, we 
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can visit at least I/7 of the pre-trimming optimal. In [T2] we filled in the gap between these 
two extremes with an 67/(5 + l)-approximation for speedup in the range 1 < s < 2 and a 47/s- 
approximation for speedup in the range 2 < s < 4. We continue to demonstrate the flexibility of 
trimming in the realm of speedup by extending these results to time windows with different lengths. 

3 The Ensemble Approach for Analyzing Performance 

Given an instance of the repairman problem on unit time windows, our previous work [12] presented 
algorithms for rational speedup s = q/r in the range 1 < s < 4 that take 0(min{r, m}r(n)) 
time, where m is the number of distinct periods. Since the approximation function is smooth and 
continuous, those algorithms work for any real speedup s, in the same range, in 0{mT(n)) time. 

The analysis of these algorithms uses a number of different service runs on trimmed time 
windows that are based on moving backward and forward along an optimal tour R*. We rely on 
averaging over a suitable ensemble of runs to establish that some run R on trimmed time windows 
does well. Because we will build on this technique and, indeed, use some of the same runs from 
our earlier work |l2j in ensembles for variable length windows, we will review our notation for 
describing these runs. 

We define unit speed to be some reference speed. Traveling with s = 1 is traveling at unit 
speed. Our results will hold whenever unit speed is no faster than the slowest speed at which 
an optimal service run is able to visit all locations during their time windows. Intuitively, this 
restriction means that we are focusing on those cases when unit speed is low enough that speeding 
up our service runs will actually give some benefit. Let R* be an optimal service run at unit speed 
originally starting at time 0. 

In our analysis we use the term racing to describe movement, forwards and backwards, along 
R* at a speedup of s times unit speed. Note that our analysis of run coverage is described on a 
period length of .5 even though, in Sect. [El we will apply this analysis to our algorithm, which uses 
three different period lengths. 

Define service run A as follows. Start run A at time t = at the location that R* has at time 
t = —0.5. Then run A follows a pattern of racing forward along R* for 1 period, racing backward 
along R* for 1 — 1/s periods, and then racing forward along R* for 1/s periods. Note that the 
pattern of movement for run A repeats every 2 periods. 

Considering the problem in which windows have length between 1 and 2, let A be an upper 
bound on the number of periods fully contained in a window. Define A^, the ^^reverse" of A, as 
follows. When run on a set of requests whose windows each fully contain at most A periods, run A^ 
starts at time t = at the location that R* has at time t = A/2. Then run A^ follows a repeating 
pattern of racing forward along R* for 1/s periods, racing backward along R* for 1 — 1/s periods, 
and then racing forward along R* for 1 period. Figure [T] shows examples of runs A and A^ with a 
speedup of 2 when A = 1. 

For the purposes of analyzing our run A, number the periods 0, 1, 2, and so on by the integer 
multiples of .5 that give their starting times. Run A repeats every 2 periods, and its coverage varies 
depending on whether the period number is even or odd. To balance this asymmetry we define 
A and A^, shifted versions of A and A^, respectively. Run A follows the same pattern as A but 
starts the pattern at time .5 at the location R* has at time 0. Run A^ follows the same pattern as 
A^ but starts the pattern at time .5 at the location R* has at time A/2 + .5. 

Our analysis will require several versions of A that have different starting points. To simplify 
notation, for any given rational speedup s = q/r, let a hop be the amount of distance traveled in 
l/(2r) time at unit speed. Let A^ be the run A moved forward A hops and let run A^ be the 
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Run A" 

I 

Optimal Run R* »~ " 
-.5 

Figure 1: Examples of runs A and with a speedup of 2 when A = 1 based on an optimal run. 
Times are labeled on the optimal run as well as runs A and A^. Segments of each run are also 
designated 5°, S^, and S'^ depending on which subset of runs they add coverage to. This subset 
naming scheme will be fully explained in Sect. HI 

run A^ moved backward A hops. Run A^ follows the same pattern of movement as run A but 
starts at time t = at the location that R* has at time t = —0.5 + A/(2r). Run A^, the reverse 
of follows the same pattern of movement as A^ but starts at t = at the location that R* 
has at t = A/2 — A/(2r). These reversed and shifted versions of A were required to establish the 
performance of our algorithms in [12J on unit length windows. Although run A with its 2-period 
repeating pattern is sufficient for those cases, we will introduce additional runs which repeat after 
3 or 4 periods in order to handle windows of longer length. 

If a service request p is serviced by a run R during the period that the time window of p has been 
trimmed into, we say that R covers p. Let S" be a subset of service requests. Define the coverage of 
5 by a run R, written cover ii{S), to be the number of requests in S covered by R divided by the 
number of requests in S. Define the coverage of S" by a set U of runs, written coveru{S), to be the 
average of coverji{S) for every run R G U. 

We will still rely on our Average Coverage Proposition from our earlier work jl2j : 

Proposition 3.1 (Average Coverage) Let {Si, S2, S3, ...Sa} be a collection of sets of service 
requests such that \J - Si gives all the service requests serviced by R* on untrimmed windows. Let 
U be a set of service runs. If 'nimi{coveri/{Si)} = fj,, then there is at least one service run R 
such that profit[R) > fi-profit{R*). 

The Average Coverage Proposition formalizes the following intuition. Let a group of service 
runs achieve some coverage over a set of requests. Let us also say that we have divided those 
requests into many different subsets, some of which may overlap, but the union of all the subsets 
is the original set of requests. If we take the subset of requests with the worst average coverage, 
some service run in the group covers a fraction of total requests no smaller than that worst average 
coverage. Otherwise, the average coverage of all subsets would be worse than the coverage of the 
worst covered subset, which is a contradiction. 

Given a way of dividing requests into subsets, we wish to prove that some set of service runs 
achieves some lower bound on average coverage. The following section will describe the algorithm 
we will use to find a service run and the analytical techniques we will use to establish a lower bound 
on its performance. This analysis will depend on carefully showing an average coverage for various 
subsets of requests defined with respect to periods induced by trimming. 
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4 Algorithm for Windows with Lengths between 1 and 2 



In [TT], we describe algorithms for s = 1 that achieve constant approximations when window sizes 
are not necessarily uniform but are close to being uniform. We extend that approach for our 
speedup problem, but specifically for windows whose lengths differ by at most a factor of 2. From 
|12j . we give an algorithm called SPEEDUP that, for unit windows, finds a run of approximately 
optimal profit at speedup s. Our approach is to modify SPEEDUP and run it with three different 
sizes of period a: .5, .75, and 1. 

For each period size, we will consider multiple starting points for a set of periods, each spaced 
.25 apart. We modify SPEEDUP appropriately so that no window starts at the beginning of a 
period, for periods of size .5, .75, or 1. This modified algorithm is called SPEEDUPW12 and is 
specified below. In this algorithm, sets of periods whose a is .5, .75, or 1 can have 2, 3, or 4 unique 
starting positions, respectively. Depending on a given period size a and starting point, a window 
will partially fill 2 periods and fully fill 0, 1, 2, or 3 periods between the 2 partial periods. For 
£ = 1,2,3 and a specified value of a, let Wi be the set of windows that completely fill exactly £ 
subintervals and partially overlap with two more of them. 



SPEEDUPW12 

PHASE 1: 

Set a to .5 and identify windows for sets Wi, W2, and W3. 
For i from to 1, 

Set the starting point for the periods to i/4. 
For j from 1 to 2, 
For k from 1 to 3, 

Trim each window in Wi to its 1** subinterval. 
Trim each window in W2 to its j*^ subinterval. 
Trim each window in VF3 to its k^^ subinterval. 
Run SPEEDUP and retain the best result so far. 

PHASE 2: 

Reset a to .75 and then identify windows for W\ and W^- 
For i from to 2, 

Set the starting point for the periods to i/4. 
For j from 1 to 2, 

Trim each window in W\ to its 1** subinterval. 
Trim each window in W2 to its subinterval. 
Run SPEEDUP and retain the best result so far. 
PHASE 3: 

Reset a to 1 and then identify windows for Wi. 
For i from to 3, 

Set the starting point for the periods to i/4. 

Trim each window in Wi to its 1** subinterval. 

Run SPEEDUP and retain the best result so far. 



When trimming, we choose from several choices of which single full subinterval to keep for each 
window. For example, for periods of length .5 and for windows in which would have three full 
subintervals, the choices for trimming will be to trim the window down to either the first, second, 
or third full subinterval. Combining these choices with the two choices associated with windows in 
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W2 and the single choice in windows in Wi yields 6 trimmings. 

The performance for speedup for windows in Wi in the range 1 < s < 4 is the same as the unit 
time window results given by our work in [12]. In Sect. O we give the performance for speedup 
for windows in W2 in the range 1 < s < 5. In Sect. [7l we give the performance for speedup for 
windows in in the range 1 < s < 6. Note that the SPEEDUP subroutine works for any real 
number 1 < s < 6; however, our analysis will assume that s is a rational number such that s = q/r. 
In the case that s is irrational, our analysis holds in the limit because the functions we find that 
bound performance in terms of s are piecewise smooth and continuous. It is worth repeating that 
our analysis uses only a period size of .5 but can still bound the performance of our algorithm with 
its three different period lengths by using careful accounting of subset coverage. 

When dealing with windows of unit length in a previous paper [12], we defined a partition of 
requests into three sets, based on which period a request was trimmed into versus which period 
an optimal run R* serviced the request in. Set T consists of requests serviced by R* in the same 
period, set E consists of requests serviced by R* in the preceding period, and set L consists of 
requests serviced by R* in the following period. For windows of length between 1 and 2, we need to 
extend this approach. We will use a superscripted S to designate that a request that was serviced 
by R* in either the first, second, third, fourth, or fifth periods with which a request overlaps. For 
requests in Wi, the sets L, T, and E will be renamed S*^, S^, and S"^, respectively. For requests in 
W2 and Ws , we will go further and use designations through and through , respectively. 
Set consists of requests serviced by R* two periods before the period into which those requests 
were trimmed, and consists of those requests serviced three periods earlier. In the same earlier 
work |12] . we further partitioned L, T, and E into Lj, Tj, and Ej for j = 1,2, . . . ,r. In a similar 
way, we will partition sets S^, S^, S"^, S^, and into r equal- length divisions, subsets Sj, Sj, 5?, 
Sj, and Sj, for any given j, j = 1, 2, . . . , r. 

Let [w, w + k) be any time window, where 1 < A: < 2. Let to be the smallest integer multiple of 
l/(2r) that is greater than w. We designate subintervals [w,uj), [uj,uj + l/(2r)), [a; + l/{2r),uj + 
2/(2r)), . . ., [uj + (4r — l)/(2r),zi; + 2) by wq, wi, W2-, ■ ■ ., Wir- For windows of length between 1 
and 2 with a given choice of period starting times, all windows fall into set Wi, W2, or W^. In 
our analysis, there is always an implied factor of 7 that accounts for the difference between the 
approximation on a tree and on a metric graph. 

We now define a procedure called CREATE- TABLE-A that describes the process of determining 
coverage for a particular speedup s for a particular run moved forward A hops. This procedure is 
a generalization of our CREATE- TABLE procedure from [12] to A > 1. (Recall that A is an upper 
bound on the number of periods fully contained in a window.) Note that CREATE- TABLE- A is not 
an algorithm that is run in the process of finding an approximation to a repairman problem with 
speedup. Rather, it provides a template that can be used to produce the tables used in analyzing 
the performance of such approximations. So that the treatment here is self-contained, we repeat 
much of our discussion of table construction from [12] , modifying it as necessary so that it can also 
handle the additional types of runs that we will introduce. 

Before CREATE- TABLE-A can be completely defined, it is necessary to explain the pattern of 
coverage generated by a run. For the kind of runs we have seen so far, type A runs, this pattern 
takes one of two forms. Let s be a rational number such that s = q/r. Type A runs repeat every 
two periods and thus can be represented with a pattern of coverage that uses a 1 to signify a subset 
covered every period and a 1/2 to signify a subset covered every other period. 

Observe the movement of type A runs, noting that, during its first period, such a run moves 
forward the same distance that an optimal run moves during q subintervals. During its second 
period of time, it moves backward the same distance than an optimal run moves during q — r 
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subintervals and then forward the same distance than an optimal run moves during r subintervals. 
Then, the pattern repeats. When s < 2, run A, during the first period in its pattern, covers q 
successive subsets as it moves forward, while in its second period covers r subsets as it moves 
forward. Note that those subsets covered as A moves backwards add nothing additional to the 
coverage. Thus, this pattern of coverage is represented as r repetitions of 1 and q — r repetitions of 
1/2. When s > 2, run A, during the first period in its pattern, covers q subsets subintervals, while 
in its second period covers q — r subsets backwards but no new subsets forward. This pattern of 
coverage is represented as q — r repetitions of 1 and r repetitions of 1/2. 



CREATE- TABLE-A(hops A) 

Let the first element of the coverage pattern be indexed at 0. 

Number the subsets through r(A + 1). 

Define function C based on the coverage pattern, such that: 



The final coverage function defined by the table is given by J^{i) + T {i). 



The values that C{i) can take on are dependent on the types of runs used. For type A runs, 
C{i) can be 0, 1/2, or 1. Runs introduced later will have a larger range of values, but it is always 
the case that < C(i) < 1. Note that the functions F and given in CREATE-TABLE-A are 
piecewise linear functions with ranges dependent on the fundamental pattern of coverage. Due to 
its construction, the combination J-{i) + J-^{i) is also a piecewise linear function and symmetrical. 
Thus, only the range < i < [r(A + l)/2j need be listed in tables. 

Although CREATE-TABLE-A gives a procedure for creating a table for a given speedup, we 
need tables expressed symbolically to prove coverage for a range of speedups. Instead of using 
specific numbers, we can leave the basic patterns of subset coverage for a given style of run (such 
as type A runs) with its shifted version in terms of q and r. By shifting this pattern r times and 
summing the results together, we account for the different alignments a time window might have 
with respect to the various subintervals. This sum is the function J^, which can be expressed as a 
piecewise linear function. Function J^^, which describes reversed runs, can be similarly described. 
To combine the two functions symbolically, we sort the end points of the subset ranges from both 
descriptions together. If, for the given range of speedups being considered, there are two end points 
which cannot be ordered, we subdivide the range of speeds so that, in each new speed range, the 
two end points in question can be ordered. Once the end points of each subset range have been 
sorted, combining the descriptions from the normal and reversed functions of the runs is achieved 
by simply summing each range. We give an example of this process in Sect. 15.11 

5 Speedup Performance for Windows in Set W2 

Recall that W2 is the set of windows that completely fill exactly two periods. We will now explore 
the speedup-performance trade-off for windows in W2 for all speedups 1 < s < 5. For set W2, our 




Define function on the same domain: 
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analysis must consider subsets wq through w^r- Throughout our analysis, we will we assign a 1 for 
full coverage and a 1/2 for half coverage of any subset. When examining the subsets for a given 
range of speedup values, the values are symmetrical around 1(^3^/2 when r is even and symmetrical 
after t(^(3r-i)/2 when r is odd. Thus, the tables and proofs we use will not list contributions for 
subset Wi where i > [3r/2j, since the contribution at Wi in these higher ranges is the same as the 
corresponding contribution at w^r-i, by symmetry. 



5.1 Speedup 1 < s < 2 for Windows in Set W2 

For the range 1 < s < 2, we can represent any rational speedup s in the form s = (r + k)/r with 



fc' 



integers r > 1 and < A; < r. For this analysis, we consider service runs A, A^, A^-k, Af!'_ 
A2r-k, and ^2r-fc' noting that A = 2. Similar to Wi for 1 < s < 2, run A covers set well, run 
A^ covers set well, and the remaining four runs plug the holes left in the spotty coverage of sets 
and 52. 

We will pass over the simpler case for A runs and use A^-k runs to give an example of how we 
construct symbolic coverage tables. For speedup s where 1 < s < 2 and A = 2, type A runs have 
a fundamental pattern of coverage of r subsets covered every period followed hy q — r = k subsets 
covered every other period. Adjusting for the offset of A = r — and making r shifted, this pattern 
yields the values for J^{i) and given in Table [2j 
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Table 2: Separate coverage functions for A^-k and in W2 when 1 < s < 2. 

Because + J-^{i) is symmetric about i = 3r/2 if r is even and after i = (3r — l)/2 if r is 
odd, we are only interested in the range < z < [3r/2j. In this range, the sub-ranges r < i <2r — k 
and 2r — k < i < L3r/2j for J- overlap with the sub-ranges r < i < r + k and r + k < i < [3r/2j 
for T^. When k < r — k, then r — k<r + k< [3r/2j. In that case, for r + k < i < [3r/2j, 
J-{i) + J-^{i) = (2r — k/2 — i) + {i — r — k/2) = r — k, as in the last interval of the middle set 
of contributions in Table [31 When k > r — k, then r — k < 2r — k < [3r/2j. In that case, for 
2r -k <i < [3r/2\, F{i) + = (r - i/2) + {i/2 - r/2) = r/2, as in the last interval of the 

middle set of contributions in Tabled! 

Similar analysis for runs A and A^ and runs ^2r-fc and produce the rest of Tables [3] and 

m The combined coverages of runs A, A^, Ar-k, ^2r-fc) ^2T--fc' °^ their respective 

shifted versions are all given in Table [3] when k <r — k and in Tabled] when k >r — k. 

Lemma 5.1 // the contributions from A and A^ are weighted by a factor of 2 and the contributions 
from Ar-k, ^^-k' ^2r-k, md are weighted by a factor of 1, the yield for all intervals is at 

least 2r + k. 
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Table 4: Contributions of runs for windows in W2 when 1 < s < 2 and k > r — k. 



Proof: We first consider the case when k < r — k, consulting Table O 

If < i < A;, then the yield for Wi is 2r + A; + i/2, which is at least 2r + k, since i > 0. 
k < i < r — k, then the yield for Wi is 2r + 3A;/2, which is greater than 2r + A;. 
r — k < i < r, then the yield for Wi is 5r/2 + A; — i/2, which is at least 2r + k, since i < r. 
If r < z < [3r/2j , then the yield for Wi is 2r + k. 
We now consider the case when k > r — k, consulting Table [H 

If < i < r — A;, then the yield for Wi is 2r + k + i/2, which is at least 2r + k, since i > 0. 
li r — k < i < k, then the yield for Wi is 5r/2 + A;/2, which is at least 2r + k, since r >k. 
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If k < i < r, then the yield for Wi is 5r /2 + k — i/2, which is at least 2r + k, since i < r. 

If r < i < 2r — k, then the yield for Wi is 2r + k. 

If 2r — k < i < [3r/2j, then the yield for Wi is also 2r + k. □ 

Theorem 5.1 For 1 < s <2, SPEEDUPW12 finds an 8j/{s + l) -approximation to the repairman 
problem on windows in set W2 in 0(min{r, m}r(n)) time. 

Proof: By Lemma [5 .11 our analysis gives no yield less than 2r + k. Since we use two copies each of 
A and and a single copy each of A^-k, A^-k^ ^2r-fc; and averaged over r different sets 

of periods, we apply the Average Coverage Proposition over 8r runs. Thus, the fraction of optimal 
profit obtained is (2r + k)/{8-fr) = ((r + A;) + r)/{8-fr) = (s + l)/(87). □ 

5.2 Speedup 2 < s < 3 for Windows in Set W2 

For the range 2 < s < 5/2, we can represent any rational speedup s in the form s = (2r + k)/r 
with integers r > 1 and < k < r — k. For this analysis, we consider service runs A, A^, ^^-2/0 
and A^_2i^, noting that A = 2. We will use three copies each of A and A^ and a single copy each of 
Ar-2k and A^_2j^. Because the generation of the tables and the case analysis needed to show the 
coverage are involved and of a similar form as Lemma [5. 1^ we have moved these details to Appendix 

Theorem 5.2 For 2 < s < 5/2, SPEEDUPW12 finds an 87/(25 — 1)- approximation to the repair- 
man problem on windows in set W2 in 0(min{r, m}r(n)) time. 

Proof: By Lemma lA. 11 the yield is at least 3r + 2k. Since three copies each of A and A^ and a 
single copy each of Ay.-2k and A^_2f, are used, averaged over r different sets of periods, the Average 
Coverage Proposition is applied over 8r runs. Thus, the fraction of optimal profit obtained is at 
least (3r + 2/c)/(87r) = ((4r + 2k) - r)/(87r) = (2s - l)/(87). □ 

Observation 5.1 Given a speedup s' > s, we can always simulate with speedup s' the runs used 
in the analysis of speedup s by introducing delays at appropriate points in each run. Thus, an 
approximation ratio of (3 at speedup s is an upper bound on the approximation ratio at speedup s' . 

By Observation l5 . 1 \ the 27-approximation for s = 5/2 implies at most a constant 27-approximation 
to the repairman problem on windows in set W2 when 5/2 < s < 3. 

5.3 Speedup 3 < s < 5 for Windows in Set W2 

For set W2 with 3 < s < 5 where s = q/r, we consider runs A, A^, and their shifts, noting that 
A = 2. When s < 4, runs A and A give full coverage for service requests in subsets of S*^ and 
and partial coverage of service requests in S"^, while runs A^ and A^ give full coverage for service 
requests in subsets of S"^ and S"^ and partial coverage of S^. When s = 4, runs A and A go further 
by also giving full coverage for service requests in subsets of S"^, while runs A^ and A^ also give 
full coverage for service requests in subsets of S^. When s > 4, runs A and A give full coverage for 
service requests in subsets of S^, S^, and 5^ and partial coverage of service requests in S^, while 
runs A^ and A^ give full coverage for service requests in subsets of S^, S"^ and and partial 
coverage of S^. Since the contributions of the A and A^ runs and their shifted versions tend to 
balance each other, we can analyze this balance between the two over all possible sets of periods 
to find a lower bound on the total profit after trimming. 
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Theorem 5.3 For 3 < s < 5, SPEEDUPW12 finds a ^j/{s — 1)- approximation for windows W2 
in 0{min{r,m}r(n)) time. 

Proof: For runs A and A, Wi earns a 1 (denoting full coverage) for each of the r sets of periods 
where < i < q — 2r, giving a total of r for each such i. For each i > q — 2r, the total decreases 
by 1/2 from the total for i — 1. For runs A^ and A^, Wi gets 1 for each of the r sets of periods 
where 3r — — 2r) < i < 3r, giving a total of r for each such i. For each i < 3r — {q — 2r), the 
total decreases by 1/2 from the total for z + 1. 

The combined contributions of runs A, A^, and their shifted versions is q/2 — r/2 for wq. 
Contributions from run A are constant, and contributions from run A^ only increase or stay 
constant for Wi where < i < r. Contributions for Wi for all runs sum to 2r for r < i < [3r/2j. 
Thus, the yield for all Wi is at least q/2 — r/2. Since two runs averaged over r different sets 
of periods arc used, the fraction of optimal profit obtained is at least {q/2 — r/2) ■ l/(2jr) = 
g/(47r) - r/(47r) = (s - l)/(47). □ 

6 New Types of Runs to Handle Windows in Set W3 

In addition to type A runs, which repeat every 2 periods, our analysis of windows in W3 defines 
type B and type C runs which repeat every 3 or every 4 periods, respectively. We define B and C 
runs only in the range 2 < ,s < 3. Let = [s] — s. 

Start run B at t = at the location that R* has at time t = —0.5. From there, run B follows a 
repeating pattern of racing forward along R* for 2 periods, racing backward along R* fov l — u/ (2s) 
periods, and racing forward along R* for v/{2s) periods. As for run C, also start it at t = at the 
location that R* has at time t = —0.5. From there, run C follows a repeating pattern of racing 
forward along R* for 2 + 2/ s periods and racing backward along R* for 2 — 2/ s periods. 

Similar to A^, we also define B^ and C^, the "reverses" of runs B and C, respectively. Both 
runs B^ and start at t = at the location that R* has at time t = \/2. Prom its starting point, 
run B^ follows a repeating pattern of racing forward along R* for v/{2s), racing backward along 
R* for 1 — I'/ {2s) periods, and racing forward along R* for 2 periods. From its starting point, run 

follows a repeating pattern of racing backward along R* for 2 — 2/s periods and forward along 
R* for 2 + 2/s periods. As with A, let Ba and Ca, respectively, be runs B and C moved forward 
A hops, and let runs B^ and C^, respectively, be runs B^ and moved backward A hops. Runs 
B/\ and Ca, respectively, follow the same patterns of movement as runs B and C but start at i = 
at the location that R* has at t = —0.5 + A/(2r). Their reverses B^ and C^, respectively, follow 
the same patterns of movement as B^ and but start at t = at the location that R* has at 
t = X/2- A/(2r). 

As there is for type A runs, there are unique patterns of coverage corresponding to type B and 
C runs and their reverses. Recall that s = q/r. Note that, for analysis of B and C runs, we choose 
the smallest values of q and r such that q + r \s even. Since the number of subsets is determined 
by r, it is necessary for B and C runs to have an even g + r in order to keep the coverage defined 
in terms of complete rather than partial subsets. A type B run moves forward during its first 
two periods of time the same distance that an optimal run moves during 2q subintervals. During 
its third period of time, it moves backward the same distance that an optimal run moves during 
q{l — v/{2s)) subintervals and then forward the same distance that an optimal run moves during 
qi'/{2s) subintervals. Then, the pattern repeats. 

We recall the list of subsets: 5?, . . . 5°, Sl . . . S^, Sj... S^, Sf . . . S^, Sf . . . Sf. When 
2 < s < 3, run B, during the first period in its pattern, covers q successive subsets as it moves 
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forward. In the second period, it covers another q subsets moving forward. Finally, in its third 
period, it covers (3g — 3r)/2 subsets backward but no new subsets forward, since the subsets covered 
forward were already covered backward. We see that only the first period in the pattern covers 
the first {q — r)/2 subsets. Then the first and third period in the pattern cover the next {q — r)/2 
subsets. All three periods cover the next r subsets. The second and third periods cover the next 
q — 2r subsets, and only the second period covers the final r subsets. This pattern of coverage is 
represented as {q — r)/2 repetitions of 1/3, (g — r)/2 repetitions of 2/3, r repetitions of 1, g — 2r 
repetitions of 2/3, and r repetitions of 1/3. 

Figure [2] gives two examples of type B runs for speedups in the range 17/7 < s < 3, the only 
range for which our analysis will employ type B runs. Observe that the run for s = 5/2 uses the 
form 10/4 in order to conform with the restriction for our analysis that q+r must be even. Unlike A 
runs which repeat every two periods, both the runs in this figure arrive at the same corresponding 
position at the beginning of every third period, namely at times 0, 1.5,3, and so on. Portions of 
runs servicing requests in S^, S^, S*^, S^, and 5^ or various subsets are identified: subsets S^, S^, 
Sl, Sf, and mapping to quarter periods of R* for s = 10/4 and subsets ^5, ^5, Sf, S^, 
and Sf mapping to a fifth of a period of R* for s = 13/5. Focusing on the example of s = 10/4 
where q = 10 and r = 4, note that, during a three-period section, subsets through 5*3 are covered 
a single time, subsets S^, Sl, and ^2 are covered twice, subsets S^, Sl, Sf, and 5^ are covered all 
three times, subsets 5| and 5| are covered twice, and subsets Sf through Sf are covered a single 
time. This pattern of 3 repetitions of 1/3, 3 repetitions of 2/3, 4 repetitions of 1, 2 repetitions 
of 2/3, and 4 repetitions of 1/3 exactly corresponds to the repeating pattern of subset coverage 
described in the previous paragraph. 




Figure 2: Examples of type B runs for two different speedups in the range 17/7 < s < 3, namely 
at 5/2 and 13/5. 

A type C run moves forward during its first two periods of time the same distance that an 
optimal run moves during 2q subintervals. During its third period of time, it moves forward the 
same distance that an optimal run moves during 2q/s subintervals and then backward the same 
distance that an optimal run moves during (1 — 2/ s)q subintervals. During its fourth period of 
time, it moves backward the same distance that an optimal run moves during q subintervals. Then, 
the pattern repeats. 

When 2 < s < 5/2, run C, during the first period in its pattern, covers q successive subsets as 
it moves forward. In its second period, it covers another q subsets moving forward. In its third 
period, it covers 2r subsets forward but no new subsets backward. Finally, in its fourth period. 
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it covers q subsets moving backward. We see that only the first period in the pattern covers the 
first r subsets. Then, the first and fourth period in the pattern cover the next q — 2r subsets. The 
first, second, and fourth periods cover the next r subsets. The second and fourth periods cover the 
next q — 2r subsets. The second, third, and fourth periods cover the next 3r — q subsets. Only 
the second and the third period cover the next q — 2r subsets, and only the third period covers the 
final r subsets. This pattern of coverage is represented as r repetitions of 1/4, q — 2r repetitions of 
1/2, r repetitions of 3/4, q — 2r repetitions of 1/2, 3r — q repetitions of 3/4, q — 2r repetitions of 
1/2, and r repetitions of 1/4. Just as with type A runs, we will use the patterns for B and C runs 
in conjunction with CREATE- TABLE-A to construct tables showing bounds on average coverage. 

Figure [3] gives two examples of type C runs for speedups in the range 2 < s < 17/7, the only 
range for which our analysis will employ type C runs. Observe that the run for s = 2 uses the form 
4/2 in order to conform with the restriction for our analysis that q + r must be even. Unlike A and 
B runs with, respectively, repeating patterns of two and three periods, both the runs in this figure 
arrive at the same position at the beginning of every fourth period, namely at times 0, 2, and so 
on. Portions of runs servicing requests in S^, S^, S'^, S"^, and or various subsets are identified: 
no separate subsets for s = 2 but subsets 5*2, 5*3, Sl, S^, 5^, Sf, Sf, Sf, and 5^ mapping to a fifth 
of a period of R* for s = 11/5. Focusing on the example of s = 11/5 where q = 11 and r = 5, note 
that, during a four-period section, subsets through 5*5 are covered a single time, subset is 
covered twice, subsets ^2 through ^5 and Sf are covered three times, subset S2 is covered twice, 
subsets 5*3 through and Sf are covered three times, subset is covered twice, and subsets 
51 through 5^ and Sf and S2 are covered a single time. This pattern of 5 repetitions of 1/4, 1 
repetition of 1/2, 5 repetitions of 3/4, 1 repetition of 1/2, 4 repetitions of 3/4, 1 repetition of 1/2, 
and 5 repetitions of 1/4 exactly corresponds to the repeating pattern of subset coverage described 
in the previous paragraph. 




Figure 3: Examples of type C runs for two different speedups in the range 2 < s < 17/7, namely 
at 2 and 11/5. 

7 Speedup Performance for Windows in Set 

We will now explore the speedup-performance trade-off for windows in W3 for all speedups 1 < s < 
6. For set W3, our analysis must consider subsets wq through W4r- As before, we will assign a 1 
for full coverage and a 1/2 for half coverage of any subset. Because of the coverage patterns of B 
and C runs, we will also assign values of 1/4, 1/3, 2/3, and 3/4 for corresponding proportions of 
coverage. For the subsets for a given range of speedup values for W3, the values are symmetrical 
around W2r- Thus, our tables and proofs will not list contributions for subset Wi where i > 2r. 
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7.1 Speedup 1 < s < 2 for Windows in Set W3 

For this analysis, we consider service runs A, A^, Aj—k, A^-ki ^2r-fc) ^2r-fe' ^3r-fci aiid A^^_f^, 
noting that A = 3. We wiU use two copies each of A and A^ and a single copy each of Ar-k, A^_j^, 
A2r-k, ^2r-fc' ^3r-fe; and We have moved the generation of the tables and the case analysis 

needed to show the coverage to Appendix iBl 

Theorem 7.1 For 1 < s < 2, our algorithm finds a 107/(5 + 1)- approximation to the repairman 
problem on windows in set W3 in 0(min{r, m}r(n)) time. 

Proof: By Lemma IB.ll our analysis gives no yield less than 2r + k. Since two copies each 
of A and A^ and a single copy each of Ar-k, ^^-k^ ^2r-fc; ^2r-fc' ^3r-k, and A^^_f^ are used, 
averaged over r different sets of periods, the Average Coverage Proposition is applied over lOr runs. 
Thus, the fraction of optimal profit obtained is at least (2r + A;)/(107r) = ((r + k) + r)/{10^r) = 
(s + l)/(107). □ 

7.2 Speedup 2 < s < 7/3 for Windows in Set W3 

For the range 2 < s < 7/3, any rational speedup s can be represented in the form s = (2r + k)/r 
with integers r > 1 and < /c < r/3. For this analysis, we consider service runs A, A^, C(3r-fc)/2) 
and C(^r-A:)/2' noting that A = 3. We will use a single copy each of A and A^ and two copies 
each of C(^3r-k)/2 and C^^_^^y2- We have moved the generation of the tables and the case analysis 
needed to show the coverage to Appendix [Cl 

Theorem 7.2 For 2 < s < 7/3, algorithm SPEEDUPW12 finds a 6'y / s-approximation to the 
repairman problem on windows in set W3 in 0(min{r, m}r(n)) time. 

Proof: By Lemma IC-H our analysis gives no yield less than 2r + k. Since 1 copy of each of A 
and A^ and 2 copies each of C(3r-fc)/2 and C^^_^y2 were used, averaged over r different sets of 
periods, the Average Coverage Proposition is applied over 6r runs. Thus, the fraction of optimal 
profit obtained is at least (2r + k)/{Q^r) = s/{'6^). □ 

7.3 Speedup 7/3 < s < 17/7 for Windows in Set W3 

For the range 7/3 < s < 17/7, any rational speedup s can be represented in the form s = {2r + k)/r 
with integers r > 1 and r/3 < k < 3r/7. For this analysis, we consider service runs A, A^, C2r-2k, 
and C^_2fc, noting that A = 3. We will use k copies each of A and ^4^ and r — k copies each of 
C2r-2k and C^_2fc. We have moved generation of the tables and the case analysis needed to show 
the coverage to Appendix [Pl 

Theorem 7.3 For 7/3 < s < 17/7, algorithm SPEEDUPW12 finds an 8-f/{s'^-4:S+7)- approximation 
to the repairman problem on windows in set W3 in 0(min{r, m}r(n)) time. 

Proof: By Lemma iD.ll our analysis gives no yield less than 3r^/4 + /c^/4. Since k copies of each 
of A and A^ and r — k copies each of C2r-2k and C^_2fc were used, averaged over r different 
sets of periods, the Average Coverage Proposition is applied over 2r^ runs. Thus, the fraction of 
optimal profit obtained is at least (3r2/4 + k'^ /4:)/{2jr'^) = ((2r + k^ - 4r(2r + k) + 7r'^) / {S-fr'^) = 
(s^ -4s + 7)/(87). □ 
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7.4 Speedup 17/7 < s < 3 for Windows in Set VI^s 

For the range 17/7 < s < 3, any rational speedup s can be represented in the form s = (2r + k)/r 
with integers r > 1 and 3r/7 < k < r. For this analysis, we consider service runs A, A^, i^^-fc+i) 
and B^_j^j^^, noting that A = 3. We will use 6r — Ak copies each of A and A^ and 3r — 3k copies 
each of Br-k+i and -B^^^^. We have moved the generation of the tables and the case analysis 
needed to show the coverage to Appendix lEl 

Theorem 7.4 For 17/7 < s < 3, algorithm SPEEDUPW12 finds a 7(23 - 7s)/(l + 3s - s^)- 
approximation to the repairman problem on windows in set W3 in 0(min{r, m}r(n)) time. 

Proof: By Lemma lE.l^ our analysis gives no yield less than 6r^ — 2rk — 2k'^. Since 6r — Ak 
copies of each of A and A^ and 3r — 3k copies each of Br-k+i and B^_i^_^^ were used, averaged 
over r different sets of periods, the Average Coverage Proposition is applied over 18r2 - Urk runs. 
Thus, the fraction of optimal profit obtained is at least (6r^ — 2rA: — 2/c^)/(7(18r^ — 14rA;)) = 
(r2 + 3r(2r + k) - (2r + A;)2)/(7r(23r - 7(2r + k))) = (1 + 3s - s2)/(7(23 - 7s)). □ 

By Observation IS . 1 1 the 27-approximation for s = 3 implies at most a constant 27-approximation 
to the repairman problem on windows in set W3 when 3 < s < 4. 

7.5 Speedup 4 < s < 6 for Windows in Set W3 

For set W3 with 4 < s < 6 where s = q/r, we consider runs A and A^, noting that A = 3. When 
s < 5, run A and its shift give full coverage in subsets of S^, S^, and S"^, and partial coverage in 
subsets of and S^, while run A^ and its shift give full coverage in subsets of 5^, S^, and S"^ 
and partial coverage in subsets of and S^. When s > 5, runs A and A give full coverage in 
subsets of S^, S^, S"^ , and and partial coverage in subsets of S^, while runs A^ and A^ give full 
coverage in subsets of S^, 5^, S"^, and and partial coverage in subsets of and . Since the 
contributions of the A and A^ runs and their shifted versions tend to balance each other, we can 
analyze this balance between the two over all possible sets of periods to find a lower bound on the 
total profit after trimming. 

Theorem 7.5 For 4 < s < 6, our algorithm finds a 47/(s — 2) -approximation for windows W3 in 
0(min{r, m}r(n)) time. 

Proof: A 1 is assigned for any subset which is covered every period, and a 1/2 is assigned for 
any subset covered every other period. For runs A and A, Wi earns a 1 for all r sets of periods 
where < i < q — 2r, giving a total of r for each such i. For each i > q — 2r, the total decreases 
by 1/2 from the total for i — 1. For runs A^ and A^, Wi gets 1 for all r sets of periods where 
4r — (g — 2r) < i < 4r, giving a total of r for each such i. For each i < 4r — (q — 2r), the total 
decreases by 1/2 from the total for i + 1. 

Now, we take the total over all sets of periods for both runs. For runs A and A, we get a total 
of r for wq. For runs A^ and A^, we get a total of r — (l/2)(4r — {q — 2r)) = q/2 — 2r. Summing 
these together, we get a yield of q/2 — r for wq. By symmetry, the total for W4r is also q/2 — r. 
Contributions from A and A are constant and contributions from A^ and A^ only increase or stay 
constant for < i < r. Contributions for Wi for all runs sum to 2r for r < i < 3r. Thus, the 
yield for all other Wi in all cases is at least q/2 — r. Since r sets of periods for the two pairs of 
runs cost a total of 2r sets of periods to average over, the fraction of profit after trimming is at 
least {q/2 — r) ■ l/(2r) = <?/(4r) — 2r/(4r) = (s — 2)/4. Multiplying the reciprocal by 7 gives a 
47/(3 — 2)-approximation. □ 
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8 Performance of SPEEDUPW12 



Now that we have characterized the performance of our speedup algorithms on windows in sets Wi, 
W2, and W3, we bound the performance of SPEEDUPW12 by combining our results as follows. 
Let R* be an optimal service run for a repairman instance with time window lengths from 1 up to 
2. Consider a new set of periods with duration .25. Partition windows into the sets H^, H^, H^, 
Hq, and H-/, such that for i = 3,4,5,6,7, a window is put in Hi if it completely contains exactly 
i of these new periods of length .25. Let the total fraction of profit in an optimal solution coming 
from windows in set Hi be hi. Thus, Yle=3 ^t- — ^• 

We use these subintervals when analyzing the performance of the algorithm run on periods 
of length .5, .75, and 1. Consider set H^ of windows, i = 3,4,5,6,7 and period length j/A, for 
j = 1, 2, 3, 4. The number of full subintervals of a window in Hi that are covered when the period 
length is j/4 is either [(^ — j + or \{i — j + depending on which set of periods is used. 

For i = 1,2,3, let fi{s)/^ be the fraction of optimal profit earned for SPEEDUPW12 applied 
to requests with windows in Wi. Recall that coverage of windows in Wi is defined using period size 
.5. To apply this coverage to the three period sizes used in the algorithm, we establish: 

Lemma 8.1 The first phase of SPEEDUPW12 yields a run R with ^■profit{R) / profit{R*) = 

Pi > /l(s)/i3 + ihflis) + i/2(s)) /J4 + f2{s)h + ikMs) + i/3(s)) he + h{s)hj . 

Proof: Windows from H^ contribute /i(s)/i3 in both sets of periods. Windows from H^ contribute 
/i(s)/i4 in one set of periods and /2(s)/i4 in the other. Windows from H^ contribute /2(s)/i5 in both 
sets of periods. Windows from H^ windows contribute /2(s)/i6 in one set of periods and f3{s)hQ in 
the other. Finally, windows from set Hg windows contribute f3{s)hj in both sets of periods. □ 

Lemma 8.2 The second phase of SPEEDUPW12 yields a run R with ^■profit{R) / profit{R*) = 

P2 > lfl{s)h3 + lfl{s)hi + /l(s)/i5 + + i/2(s)) he + + |/2(S)) hT. 

Proof: Windows from H^ contribute |/i(s)/i3 in one set of periods and nothing in the other two. 
Windows from Hq contribute \fi{s)hi in two sets of periods and nothing in the other one. Window 
from Hj contribute \fi{s)h^ in all three sets of periods. Windows from H^ contribute |/i(s)/i6 in 
two sets of periods and ^/2(s)/i6 in the other one. Windows from Hg contribute g/i(s)/i7 in one 
set of periods and ^f2{s)h-j in the other two. □ 

Lemma 8.3 The third phase of SPEEDUPW12 yields a run R with ^■profit{R) / profit{R*) = 
P3 > \fi{s)hi + |/i(s)/i5 + |/i(s)/i6 + fi{s)hT. 

Proof: Windows from H^ contribute nothing in all four sets of periods. Windows from H^ con- 
tribute \fi{s)hi in one set of periods and nothing in the other three. Windows from H^ contribute 
|/i(s)/i5 in two sets of periods and nothing in the other two. Windows from ifg contribute |/i(s)/i6 
in three sets of periods and nothing in the other one. Windows from Hg contribute \fi{s)h'j in all 
four sets of periods. □ 

From Lemmas 18. H 18.21 and 18.31 we isolate the coefficients bi of the variables hi, for £ = 3, 4, 
5, 6, 7. Weighting them by x, y, and z to correspond to those lemmas, respectively, leads to the 
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following definitions of five functions of x, y, z, and s. 
h = h{s)x + ^fi{s)y 

h = (^lfi{s) + lf2{s)y+pi{s)y + ^Ms)z 

1 

h = f2{s)x + fi{s)y + -fi{s)z 

h = Q/2(^) + ^/3(^))x+Q/i(.) + ^/2(.))y + ^/i(.)z 
br = his)x+(^^his) + ^his)^y + fi{s)z 

Theorem 8.1 In 0(min{r, m}r(n)) time, SPEEDUPW 12 finds a run R such that ^■profit[R) 
/{profit{R*)) > max{ p \ p < be for £ = 3,4,^,^,7, x + y + z < 1, x > 0, y > 0, z > }. 

Proof: By Lemmas 18.11 18.21 and 18.31 profit{R)/{profit{R*)j) > max{/>i, p2, pa}. Then, for any 
convex combination of pi, p2, ps, (i.e., x,y, z > and x + y + z = 1), we have 

profit (R) \ ^ A 

> max < mm {pix + p2y + psz} } 



profit{R*)j x+y+z=l [^£{3,4,5,6,7} and hi=l 
x,y,z>0 

Thus, the expression for be is given by setting he = 1 and hi = 0, where i ^ i, and summing 
pi, p2, and ps, weighted by x, y, and z, respectively. In this way, we can account for a problem 
instance being dominated by any set He for £ = 3, 4, 5, 6, 7. No problem instance will be worse than 
a convex combination of all the bounds. 

Algorithm SPEEDUPW12 runs SPEEDUP a total of 12 times in the first phase, 6 times in the 
second phase, and 4 times in the third phase, for a total of 22 times. We have shown that the 
running time of SPEEDUP is 0(min{r, m}r(n)). □ 

We give a description of /i(s), f2is), and fais) in Tabled Function /i(s) comes from our work 
in [12]. Function /2(s) comes from Sect. [5l Function fsis) comes from Sect. [71 

We produced the results in Tabled] by solving the linear programs of Theorem 18. II for particular 
values of s within each range, inferring the pattern for each range, and then proving the inferred 
pattern. Note that all but one of the reciprocals of the resulting ratios in terms of s are nonlinear 
functions! 

Theorem 8.2 For speedup s in the range 1 < s < 6 and window lengths between 1 and 2, algo- 
rithm SPEEDUPW12 produces a service run R for the repairman problem with approximation ratio 
profit{R*) / profit{R) upper-bounded as in Table\^ 

Proof: For each possible speedup range, we show that 7 times the convex combinations of the 
functions given in Table [5] are never less than the reciprocals of the approximation ratios listed in 
Tabled] 

50 6 , 17 , , , 26s + 26 
When 1 < s < 2, choose x = — , y = — , and z = — . Ihen, 03 = 64 = . . . = 67 = . 

wv, V. + 3s -3s^ + 9s 4s^ - 6s + 3 

When 2 < s < 7/3, choose x = — ^ , y = — , and z = — ^ ■ 

- - ' ' 7s2 + 6s + 3 ' ^ 7s2 + 6s + 3 ' 7s2 + 6s + 3 
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his) > 



f2{s)> { 



his) > 




( is + l)/8 
(2s - l)/8 
1/2 
I is - l)/4 

(s + l)/10 

s/6 

(s2 -4s + 7)/8 

(1 + 3s - s2)/(23 - 7s) 

1/2 

is - 2)/4 



1 < 

2 < 



1 < s 

2 < s 
1 < s _ 

3 < s < 



1 

2 

7 

3 
17 
7 

3 



< 2 

< 4 

< 2 



< 
< 
< 
< 
< 
< 



< 5 

— 2 

< 3 



< 
< 
< 
< 
< 
< 



7 
3 

iZ 

7 
3 

4 

6 



Table 5: Lower bounds on fractions of optimal profit collected for the sets Wi, W2, and W3, ignoring 
the factor of 7. 



5s3 + 6s2 

Then, 63 = . . . = 67 = 



28s2 + 24s + 12 ■ 



, 4s2 + 2s 3s3-18s2 + 27s 

When 7/3 < s < 17/7, choose x = — ^ — — ^ — -, y 



-53 + 1052-38 + 2"" -53 + 1052-35 + 2' 

4s3 - 24s2 + 32s - 2 , , s^ - 25^ + lls^ 

and z = ^ . Then, 63 



-53 + 1052 -3s + 2' ' -453 + 4052-125 + 8' 

..ru 1 ^ / ^ . /n r, 1453 - 39s2 - 23s -953 + 54^2 _ gj^ 

When 17/7 < s < 5/2, choose x = — — — — — , y 



1753 - 4352 - 355 - 23 ' ^ 1753 - 4352 - 355 - 23 ' 

1253-5852 + 695-23 , , lls^ - 2I53 - 5O52 

and z = „ — — ^5 — — . Then, 63 = . . . = 67 



1753 - 4352 - 355 - 23 ■ ' ■ ■ ■ 6853 - 17252 - 140s - 92 ' 

28s3 - 120s2 + 92s 33s3 - I8952 + 2645 

When 5/2 < s < 3, choose x = ——5 , „ — — — — , y 



7353 - 40952 + 6685 - 368 ' ^ 7353 - 409s2 + 6685 - 368 ' 

1253 - 100s2 + 312s - 368 , , 395"^ - 183s3 + 180s2 

and z = — — — — ■. Ihen, 63 = . . . = 07 = 



7353 -409s2 + 6685 -368' ' "' 29253 - 1636s2 + 26725 - 1472 ' 

..rr. 0/ u 2s2 + 2 -3s2 + 12s , 4s2 - 12s + 4 
When 3 < 5 < 4, choose x = — -, y = --^ — -, and z 



Then, 63 = . . . = 67 = 



352 + 25 + 4"" 3s2 + 2s + 4' 3s2 + 2s + 4 

53 + 6s2 



1252 + 85 + 16' 



When 4 < 5 < 5, choose x = — y = — 5f_lLi^ g^^j z = 0. 
- - ' -s + 16' -5 + 16 ' 
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rr.! , , s + 4 , , , - 6s + 45 s + 4 

Then, bs = br = — — , and 64 = Oe = — — — > — -, whenever s < 16. 

—s + 16 —4s + 64 —s + 16 

s*^ — 6s + 45 s + 4 

This follows since whenever s < 16, > holds if and only if (s — 5)^ + 4 >0, 

—4s + 64 — s + 16 

which is always true. 

1^- 11 u Au s^-8s + 37 s + 4 , 

J^mally, bound 05 = > whenever s < 16. 

—2s + 32 — s + 16 

s^ — 8s + 37 s + 4 

This follows since whenever s < 16, > holds if and only if (s — 5)^ + 20 > 0, 

—2s + 32 — s + 16 

which is always true. 

When 5 < s < 6, choose x = , y = — — , and z = 0. 

- - -3s + 26'^ -3s + 26 

s - 14 , , , 2s - 20 s - 14 , 

Then, o-i = bv = , and 04 = or = > whenever s < 6. 

3s - 26 3s -26 -3s -26 

s — 14 

Finally, 65 = 1 which is at least whenever s < 6. □ 

3s — 26 



9 Conclusion 

This paper has demonstrated the surprising versatility of the technique of trimming. Even with time 
windows whose lengths are not all the same, it is possible to simplify the structure of many time- 
constrained route-planning problems and apply an ordering that allows dynamic programming to 
work well. For unrooted problems, the cost of this additional order is at most a constant reduction 
in the profit a run can earn. We have extended results from our previous paper [12] so that we can 
characterize the way in which this reduction in profit can be offset, in part or in whole, by speedup 
over a hypothetical optimal benchmark when the lengths of time windows are not all uniform. The 
key idea needed for this extension is to consider a diverse set of trials with a number of different 
period lengths for trimming and then choose the best result among all those found. This approach 
makes trimming adapt to various distributions of window lengths. 

We have given techniques that achieve an approximation ratio parameterized only by speedup 
when the ratio between the longest time window and the shortest time window is no greater than 
2, but these techniques can be extended to other ranges of time window lengths. For the general 
case, in which the ratio between the longest and the shortest time windows is D, the approximation 
ratio will worsen by a factor of log2 D, using an approach similar to the one we used in [11] for 
general length time windows without speedup. 

It is worth mentioning that we have achieved approximation bounds for a few specific ranges 
of s which are slightly better than the ones listed in Table [TJ While trying to accommodate these 
ranges into a coherent scheme, our analysis became so much more complex that we chose to give 
a more complete and readable presentation of results which are nearly as strong as the best we 
found. The fact that better values are possible shows that there is potential in these techniques. 
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A Coverage of Windows in Set W2 when 2 < s < 3 

Analysis of set W2 when 2 < s < 3 is done by considering service runs A, A^, Ar-2k, and Af!'_2f^, 
noting that A = 2. The combined coverages of runs A, A^, Ay._2k, ^^_2fci ^^^'^ of their respective 
shifted versions are given in Table[6]when k < r — 2k and in Table[7]when k > r — 2k. The combined 
coverage of the pair A and A^ is exactly the same for all values of k and are only listed in Table [6j 
These and all other tables in the Appendices are generated by using CREATE- TABLE-A for each 
different run type, using the appropriate value of A (2 for W2 or 3 for W3), and specific values for 
A determined by the number of hops each run has been moved. 

r 0< i <k 

r + ^k — k < i < r — k 

^r + k r-k< i <[^\ 



Combined contributions 
for A and A^ 



Combined contributions 
for Aj.^2k and A? 



2k 



( 2k + i 

^k + |i 
r + ^i- 
k 



\k 



< 
k < 
r-2k< 
r — k < 



< k 

<r-2k 

< r — k 

< [^1 



V 2' ' — " — L 2 J 

Table 6: Contributions of runs for windows in W2 when 2 < s < 5/2 and k <r 



2k. 



Combined contributions 



( 2k + i 0< i <r-2k 

r r — 2k < i < k 

k k < i < r 

^2' ^ r — k < i < \_ 

Table 7: Contributions of runs for windows in W2 when 2 < s < 5/2 and k > r 



for Ar-2k and A, 



R 

r-2k 



k 



3r 



2k. 



Lemma A.l // the contributions from A and A are weighted by a factor of 3 and the contributions 
from Ar-2k o-^d A^2fe '^^^ weighted by a factor of 1, the yield for all intervals is at least 3r + 2k. 

Proof: We first consider the case when k <r — 2k, consulting Table O 

If < i < A:, then the yield for Wi is 3r + 2/c + i, which is at least 3r + 2k, since i > 0. 

k < i < r — 2k, then the yield for Wi is 3r + 3k, which is greater than 3r + 2k. 
li r — 2k < i < r — k, then the yield for Wi \s Ar + k — i, which is at least 3r + 2k, since i < r — k. 
li r — k < i < [3r/2j , then the yield for Wi is 3r + 2k. 

We now consider the case when k > r — 2k, consulting Tables [6] and [71 The algebra for the 
cases when 1 < i < r — 2k, k < i < r — k, and r — k < i < L3r/2j gives exactly the same results as 
the first, third, and fourth ranges from the previous part of the proof. r — 2k < i < k, then the 
yield for Wi is 4r > 3r + 2k, since r > 2k when 2 < s < 5/2. □ 
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B Coverage of Windows in Set W3 when 1 < s < 2 

Analysis of set W3 when 1 < s < 2 is done by considering service runs A, A^, Ar-k, ^^^-fc, ^2r 



-fc, 



A2r-k, A^r-k^ ^3r-fc, AS--*.' ^nd all of 



^3r-fc, and A§j._f^, noting that A = 3. 
The combined coverages of runs A, A^, Ar-k, ^r—k^ ^-zr-k, ^2r-fc' ^ir-k, ^zr-k 
their respective shifted versions are given in Table [8] when k < r — k and in Table [9] when k > r — k. 
The combined coverages of the pair A and A^ and the pair Aj.-k and A^_i^ are exactly the same 
for all values of k and are only listed in Table \8\ 



Combined contributions 
for A and A^ 



r — 
r + 



0< i <k 

k < i < r 

r < i <r + k 

r + k < i < 2r 



Combined contributions 
for Ar-k and A^!'_, 



Combined contributions 
for A2r-k and A^^_i^ 
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for A^r-k and ^f^.^ 





k + i 
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< r + k 




2z - 2r + ifc 
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+ k 
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i 


<2r-k 




+ r — k 


2r 


- k 
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i 


< 2r 
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i 


< k 




i — 2^ 
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i 


< r 
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2^+2'^" 2^ 
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< r -\- k 




2r + /c — z 
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+ k 
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i 


<2r-k 




2A; 


2r 


- k 


< 


i 


< 2r 



Table 8: Contributions of runs for windows in W3 when 1 < s < 2 and k < r — k. 



Lemma B.l /f the contributions from A and A^ are weighted by a factor of 2 and the contributions 
from Ar-k, ^^-k> ^2r-fc; ^2r-k' ^3r-A:> c-^d are weighted by a factor of 1, the yield for all 

intervals is at least 2r + k. 

Proof: We first consider the case when k < r — k, consulting Table [8l The algebra for the cases 
when < i < r gives the same results the first, second, and third cases in Lemma I5.H at least 
2r + fe in each case. If r < i < 2r, then the yield for Wi is 2r + k. 

We now consider the case when k > r — k, consulting Tables [8] and [H The algebra for the cases 
when < i < r gives the same results as the proof of Lemma |5. II for k > r — k, at least 2r + in 
each case. If r < i < 2r, then the yield for Wi is again 2r + k. □ 
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Combined contributions 
for A2r-k and A^^_f, 



Combined contributions 
for A^r-k and Af^_f^ 



r + k 



i + r — k 



Ik 



2^ + ^ - - 



2A; 



2'" 



< 

r — k < 
r < 
2r - /c < 
r + k< 

< 

< 
r < 
2r-k < 
r + k < 



< r — k 

< r 

< 2r - /c 

< r + fe 

< 2r 

< A; 

< r 

<2r-k 

< r + k 

< 2r 



Table 9: Contributions of runs for windows in W'i when 1 < s < 2 and k > r — k. 



C Coverage of Windows in Set W3 when 2 < s < 7/3 

Analysis of set W3, when 2 < s < 7/3 is done by considering service runs A, A^, C(3r_fc)/2) and 
C(^,._fc)/2' looting that A = 3. The combined coverages for these runs are listed in Table [TO] assuming 
that r+k is even. When r+k is not even, we can achieve an identical speed by multiplying both by 2. 



Combined contributions 
for A and A^ 



r + i/c — |i 



2r 
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k < 
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i < k 
i < 2r 
i < 2r 
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for C^3r-k)/2 and Cg,,„;.)/2 



+ ^k + 
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i(r -3A:) < 
i(r- A;) < 
(3r - 3k) < 
i(3r - A;) < 



<l(r-3A:) 
<i(r-A;) 
< i(3r-3A;) 
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i(3r - A:) 
< 2r 



Table 10: Contributions of runs for windows in W3 when 2 < s < 7/3. 



Lemma C.l If the contributions from A and A^ are weighted by a factor of 1 and the contributions 
from C(3r_fc)/2 (^^d C^^_i^y2 are weighted by a factor of 2, the yield for all intervals is at least 2r+k. 

Proof: Consulting Table \W\ we first consider the case when 5k < r, which implies k < {r — 3k) /2. 

If < i < k, then the yield for Wi is 2r + Ak + i, which is at least 2r + A;, since i > 0. 

If A; < i < (r — 3k)/2, then the yield for Wi is 2r + 3A;/2 + i/2, which is greater than 2r + k. 

If (r — 3A;)/2 < i < {r — k)/2, then the yield for Wi is 5r/2 — i/2, which is at least 2r + k, since 
i < (r — k)/2 and r > 3A;. 

If (r — A;)/2 < z < (3r — 3k) /2, then the yield for Wi is 9r/4 + k/4, which is at least 2r + k, since 
r > 3A;. 
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If (3r — 3k)/2 < i < (3r — k)/2, then the yield for Wi is 3r/2 + k + i/2, which is at least 2r + k, 
since i > (3r — 3k) /2 and r > 3/::. 

If (3r — k) /2 < i < 2r — k, then the yield for is 3r + k/2 — i/2, which is at least 2r + k, since 
i<2r -k. 

If 2r — k < i < 2r, then the yield for Wi is 2r + A;. 
Next, we consider the case when 5k >r, again consulting Table [TOl 
The first case gives the same result as above but for the range < i < (r — 3k) /2. 
If (r — 3k) /2 < i < k, then the yield for Wi is 5r/2 — k/2, which is at least 2r + k, since r > 3k. 
The third case gives the same result as above but for the range k <i <{r — k)/2. The fourth, 
fifth, and sixth cases above give identical results when 5k > r. □ 



D Coverage of Windows in Set W3 when 7/3 < s < 17/7 

Analysis of set W3 when 7/3 < s < 17/7 is done by considering service runs A, A^, C2r-2k^ and 
^2r-2fc' noting that A = 3. The combined coverages for runs A and A^ are listed in Table [TOl and 
the combined coverages for runs C2r-2k and C|^_2^, are listed in Table [TT| assuming that r+k is even. 



Combined contributions 
for C2r-2k and C^_2fc 



\r + \k + \i 
jr + -gk + 2^ 



Ir + k + li 



< 
r-2k < 

r < 
r + k < 
2r -k < 
r + 2k< 



<r-2k 

< r 

< r + A; 
<2r-k 
<r + 2k 

< 2r 



Table 11: Contributions of runs for windows in when 7/3 < s < 17/7. 



Lemma D.l If the contributions from A and A^ are weighted by a factor ofk and the contributions 
from C2r-2k o.nd C2r-2k ^'"^ weighted by a factor of r — k, the yield for all intervals is at least 
3rV4 + A;V4. 

Proof: Consult Tables [inland ni 

If < i < r - 2A:, then the yield for wt is 3r^/4 + k'^/4. 

li r — 2k < i < k, then the yield for Wi is r^/2 + 3rA;/4 + (r — k)i/4:, which is greater than 
3r2/4 + k'^/i, since i>r - 2k. 

If k < i < r, then the yield for Wi is r^/2 + 3rA;/4 + k'^ /A + (r — 3k)i/4:, which is at least 
3r^/4 + A;^/4, since i < r and r < 3k. 

li r < i < r + k, then the yield for Wi is r^/4 + rA; + /c^/4 + (r — 2k)i/2, which is at least 
3r^/4 + A;^/4, since i > r. 

li r + k < i < 2r — k, then the yield for Wi is 3r-^/4 + rk — A;^/4 — ki/2, which is at least 
3r^/4 + A;^/4, since i <2r - k. 

If2r — k<i<r + 2k, then the yield for Wi is r^/4 + 3rA;/4 + (r — k)i/A, which is at least 
3r^/4 + A;^/4, since i>2r - k. 

li r -\- 2k < i < 2r, then the yield for Wi is r'^/2 + rk — /c^/2, which is at least 3r^/4 + A;^/4, 
since {r"^ /2 + rk - k^ /2) - {Sr'^/A + k'^/A) = ((r - A:)/2)((3A; - r)/2) > when A; < r < 3A;. □ 
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E Coverage of Windows in Set W3 when 17/7<s<3 

Analysis of set W3 when 17/7 < s < 3 is done by considering service runs A, A^, Br-k+i, and 
noting that A = 3. The combined coverages for runs A and A^ are hsted in Table \T0\ and 
the combined coverages for runs B^-k+i and B^^j^^^ are listed in Table [T2] when 17/7 < s < 5/2 
and in Table [T3l when 5/2 < s < 3, assuming in both cases that r + A: is even. 



Combined contributions 
for Br-k+i and -B^_fc+i 



|A; + I? 

— ir + k + i 

2 



3^^+ 3^ 



< 

r — k < 

i(r + A:) < 

i(3r-3A;) < 

2r - A; < 

i(3r + /c) < 



< r — k 

< ^{r + k) 

< i(3r - 3A:) 

< 2r - A; 
<i(3r + A;) 

< 2r 



Table 12: Contributions of runs for windows in W3 when 17/7 < s < 5/2. 



Combined contributions 

for Br-k+i and B^_j^j^-^ 



Ik + li 
-ir + A; + i 

iA; + z 

Ir + ^k 



< 

r — k < 
i(3r - 3A:) < 
i(r + A;) < 
2r - A: < 
i(3r + A;) < 



< r — k 

< i(3r -3A;) 
<i(r + fc) 

< 2r - A; 

< ^(3r + A:) 

< 2r 



Table 13: Contributions of runs for windows in VF3 when 5/2 < s < 3. 



Lemma E.l /f the contributions from A and A^ are weighted by a factor of 6r — 4A; and the 
contributions from B^^k+i o.nd B^_i_^^ are weighted by a factor of 3r — 3k, the yield for all intervals 
is at least 6r^ — 2rA; — 2A;^. 

Proof: In the case that 17/7 < s < 5/2, consult Tables HO] and [HI 

If < i < A;, then the yield for Wi is 6r^ — 2rA; — 2A;^ + 2(r — k)i, which is at least 6r^ — 2rk — 2/c^, 
since r > k. 

If k < i < r — k, then the yield for Wi is 6r^ -j-rk — 4A:^ — ir = 6r^ — 2rA: — 2A:^ + {r — k — i)r + 
(r — 2A:)A; + (3A: — r)r, which is greater than 6r^ — 2rA; — 2A:^, since i < r — k, r > 2k, and k > r/3. 

Ifr — A;<i< (r + A;)/2, then the yield for Wi is 5r^ + 3rA; — 5A;^ — ik = 6r^ — 2rA; — 2A;^ + ((r + 
A;)/2 — i)k + (r — 2A;)7A;/4 + (7A; — 3r)r/3 + 5rA;/12, which is greater than 6r^ — 2rA; — 2A:^, since 
i<{r + k)/2,r> 2k, and k > 3r/7. 

If (r + A;)/2 < f < (3r - 3A;)/2, then the yield for Wi is 9r'^/2 + 3rA: - 9A;2/2 + (r - 2k)i = 

- 2rk - 2k'^ + (i - (r + k)/2){r - 2k) + (r - 2k)7k/i + (7A; - 3r)r/3 + 5rA;/12, which is greater 
than 6r^ — 2rA; — 2A;^, since i > (r + A;)/2, r > 2A;, and k > 3r/7. 

If (3r — 3A;)/2 < i < 2r — k, then the yield for Wi is Gr"^ — 3k'^ — ik, which is at least 6r'^ — 2rk — 2k^ , 
since i < 2r — k. 
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If 2r — k < i < (3r + k) /2, then the yield for Wi is 2r^ + Ark — Ak"^ + 2(r — k)i^ which is at least 
6r^ — 2rk — 2k'^, since i > 2r — k. 

If {3r + k)/2 <i <2r, then the yield for Wi is 5r'^ + 2rk-5k'^ = dr"^ -2rk -2k'^ + {r -2k)3k/2 + 
{7k — 3r)r/3 + rk/Q, which is at least 6r^ — 2rk — 2/c^, since r >2k and k > 3r/7. 

In the case that 5/2 < s < 3, consult Tables [TOl and [T3l For this range (3r — 3k)/2 < (r + A;)/2. 

If r — A; < i < (3r — 3A;)/2, then the yield for Wi is 5r^ + 3rA; — 5k'^ — ik = 6r^ — 2rk — 2k^ + 
((3r - 3A;)/2 - i)A; + (r - /c)3A;/2 + (2A; -r)r, which is at least Gr"^ -2rk -2k'^, since i < (3r-3k)/2 
and k < r < 2k. 

If (3r — 3k) /2 < i < (r + /i:)/2, then the yield for Wi is 13r^/2 — 7k'^/2 — ir, which is at least 
6r^ — 2rk — 2A;^, since i < {r + k)/2 and r > k. 

If {r-\-k)/2 < i < 2r — k, then the yield for Wi is 6r^ — 3k'^ — ik, which is at least 6r^ — 2rk — 2k'^, 
since i < 2r — k. 

All other ranges are identical to some yield when 17/7<s<5/2. □ 
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